Long DNA molecule as a pseudoscalar liquid crystal 
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We show that a long DNA molecule can form a novel condensed phase of matter, the pseudoscalar 
liquid crystal, that consists of aperiodically ordered DNA fragments in right-handed B and left- 
handed Z forms. We discuss the possibility of transformation of B-DNA into Z-DNA and vice 
versa via first-order phase transitions as well as transformations from the phase with zero total 
chirality into pure B- or Z-DNA samples through second-order phase transitions. The presented 
minimalistic phenomenological model describes the pseudoscalar liquid crystal phase of DNA and the 
phase transition phenomena. We point out to a possibility that a pseudoscalar liquid nano-crystal 
can be assembled via DNA-programming. 
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DNA is a biopolymer that provides with the basic mech- 
anisms of biological information processing. Schrodinger 
was the first who recognized that a gene or the whole 
chromosome thread is an aperiodic crystal fl| . It consists 
of two types of base pairs (AT and CG nucleotide pairs) 
which are placed in a random order along the molecular 
chain. Meanwhile natural DNA exists as a double helix 
bound state [^]. Each base in one strand is connected to 
a base in the other strand by hydrogen bond link. Suc- 
cessive bonds are slightly angled to each other leading to 
the famous double-helix form of DNA. Depending on its 
environment and mechanical strain, the DNA double he- 
lix may have several different helical forms, of which the 
most common are right-handed A- and i?-DNA 0], and 
the most peculiar is the left-handed Z-DNA Q. Several 
different helical conformations of DNA may even coexist 
as domains of the same long molecule Q as well as trans- 
form into each other. 

Previously it was pointed out Q that different con- 
densed matter phases may be formed inside long poly- 
mers, including a liquid crystal phase. In the present pa- 
per we will concentrate on the formation of a novel liquid 
crystal phase within a long DNA molecule. The main 
idea here is that the pieces of the molecule may have dif- 
ferent helicity (being in B or Z forms). This situation 
is similar to the case of the pseudoscalar liquid crystal 
that was hypothetically introduced by Zel'dovich back in 
1974 0. The peculiarity of such a phase is that a lo- 
cal right-handed helix may transform into a left-handed 
one and vice versa meanwhile interacting with other he- 
lical parts of the molecule. Change of the temperature 
leads to an ordering that happens via a spontaneously bro- 
ken chiral symmetry phase transition. The ordered phase 
thus represents the pseudoscalar liquid crystal The 
order parameter for liquid crystals is the orientation of its 
molecules and for a pseudoscalar liquid crystal the order 
parameter is the chirality ^ . In our case it is the local dif- 
ference between the number of oppositely-handed helices. 
A substantial difference from the other liquid crystals Q 
consisting of chiral molecules is that the helices in pseu- 
doscalar crystals can change their chirality under external 
conditions. 

Consider a DNA molecule as a sequence of pieces con- 
taining right- and left-handed helices, being in one of the 



two B OT Z forms with opposite chiralities. The pieces 
interact with each other that may lead, as we show be- 
low, to phase transitions bringing the molecule, e.g., to 
a pure left-handed Z-DNA. The proposed mechanism can 
explain how one of the structures, say the Z-form, emerges 
during intracellular processes. For the quantitative pur- 
poses we will employ a combination of Ising-type nearest- 
neighbor interactions and random interactions of distant 
along the chain parts which may come close. The long- 
range interactions may occur due to random bending of 
the DNA molecule caused, e.g., by looping proteins \w^ . 

We now present the following simple phenomenological 
model for description of pseudoscalar liquid crystal phase 
of DNA. This is a statistical mechanical model defined on 
a one-dimensional lattice with each site corresponding to 
a rung of the ladder. A spin variable Xi is associated with 
each site i, where we choose Xi = ~1 to correspond to the 
energetically more favorable right-handed helical rotation 
(that means the particular base pair belongs to the B-form 
segment of the molecule) and Xi ~ +1 for the base pair 
in the left-handed Z-form. Notice that one-dimensional 
Ising models with chiral variables were applied to s tudy 
random copolymers and bent-core liquid crystals in [llj. 
Their study can be easily extended to account for long- 
range interactions, e.g., via loop formation in the case of 
copolymers. 

It is reasonable to assume that the neighbor base pairs 
will try to get into the same form and thus we will take 
the nearest-neighbor interactions to be ferromagnetic. Be- 
sides the nearest-neighbor interactions we will take into 
account random interactions between base pairs which are 
distant along the chain but come close in space. We as- 
sume an arbitrary folding of the DNA molecule so that 
any two base pairs may get connected, e.g., via the loop- 
ing proteins [iO|. The proposed model has the following 
Hamiltonian 

N-l N N N 

H = XiXi+i - 7 X! "^ioXiX] - h^Xi + a^Jij 

i—1 i<j i—1 i<j 

(1) 

where g > is the coupling parameter of nearest-neighbor 
interactions responsible for the torsional stiffness of the 
polymer; Jy are the link variables, taking values and 
1 when the i and j nodes are uncoupled or coupled, e.g., 



by looping proteins, correspondingly; 7 > is the en- 
ergy of interaction between two distant base pairs cou- 
pled via appeared link; h is the half of the energy differ- 
ence between the states with opposite helicity; a is the 
energy of formation of a link connecting i and j sites. We 
introduce and will use another parameter c defined via 
c/{N — c) = e~"^, where /3 is the inverse temperature. 
The ratio can be roughly treated as the probability of 
a link formation (see [12i] for more rigorous formulations 
and details of a related model that describes a network 
of fluctuating links). We will assume sparse connectivity 

^ 1 as we suppose that the number of uncoupled dis- 
tant base pairs that tend to get coupled is much less than 
the number of base pairs N. 

The proposed model can be reduced to an effective 
mean-field Hamiltonian with the link variables being elim- 
inated. Indeed, we proceed with calculation of the par- 
tition function Z = Tr^^Trj exp [— where Tr^ = 

n!Ii Ex.=±i and Trj = J],^ Ej.,=o,i are traces over the 
corresponding variables. Taking the trace Trj we obtain 
the following expression for the partition function 

Z = Tr;^e'^^^^^'i''^''^'+'+''S^=i^'] 

N 

x[](l + e-^"+'3^x.x.) 



Recalling the definition of c and that c ^ iV, we get 
for the last product factor exp J2^j e^''^^^^ . Not- 
ing then the identity (r — ±1) e^'^^ = 60 + biT, br = 
taiih^ (Pj) cosh(/37) (r — 0, 1) we can now rewrite the par- 
tition function as 



exp[-^sinh/37^XiXj] 



N 



Thus we have arrived at the partition function with the 
following effective Hamiltonian 



iV-l 



N 



N 



H 



(2) 
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where the coupling parameter is 7' — ■^sinh/37. The 
first term is responsible for interactions between nearest- 
neighbor base pairs. It assures that neighbor pairs tend 
to get the same value of chirality. The second term in the 
Hamiltonian describes the effective mean-field interaction 
between base pairs. The last term is the half of the energy 
difference between the right- and left-handed helices at the 
sites of the chain. 



We now introduce the order parameter = J2iL 

1 



1 Xi 

the 



which takes the values ^ = — 1 and ^ — 1 for 
pure samples of B and Z forms of DNA, respectively. 
To calculate the partition function for the above mean- 
field Hamiltonian we use the relationship X]i<j XiXj — 
^{J2iLi Xi)'^~^^i the Hubbard-Stratonovich transforma- 
tion e^"(^.=i^') 
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FIG. 1: Chirality vs temperature. The parameters are g = 0.2 
kcal/mol, 7 — 0.02 kcal/mol, c = 10 and h = 0.01 kcal/mol. 
There is a second-order phase transition at the critical tem- 
perature Tc ~ 0.6 kcal/mol. 



the expression for the partition function of the one- 
dimensional (ID) Ising model ^Ts^. Then the partition 
function takes the form Z cx dfj.e^^^^'-^'' with the 
effective free energy /(/i) given by 



^sinh^/3(/i + &^) + e-4/3s] 



(3) 



where we have defined b — csinh/37. 

To further analyze equilibrium properties of the sub- 
stance under consideration the extreme values of fi are to 
be obtained via /'(/i) = that leads to the following state 
equation 



M = 



sinh (3{h + 5/i) 



Y/sinh2/3(/i + 6/i) + e-4/39 



(4) 



The dependence of the order parameter fi on temperature 
T is presented in Fig[T] For the zero-field case {h — 0) one 
would have a second-order phase transition. The molecule 
spontaneously transforms into a structure that has non- 
zero total chirality. It chooses between B and Z forms. 
The critical temperature of the phase transition is deter- 
mined via (3cb ■ e^^"^ = 1. If the field is non-zero then the 
system goes through the phase transition to get the form 
that is favored by the sign of the field. The field h de- 
pends on the external parameters that are determined by 
environmental conditions such as pH value, salt concen- 
tration, and other chemical as well as mechanical factors 
(e.g., locally applied forces and torques). 

For a fixed constant temperature there can be a first- 
order phase transition for chirality versus the field as 
demonstrated in FigO The phase transition transforms 
the right-handed B-DNA into left-handed Z-DNA at the 
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FIG. 2: Chirality vs the base pair binding energy for the pa- 
rameters g = 8.5 kcal/mol, 7 — 0.02 kcal/mol, c — 10, and 
fcsT — 0.6 kcal/mol. The first-order B — Z phase transi- 
tion occurs at the critical value he — 0. Dashes indicate the 
metastable states. 



critical point he — where the external parameter changes 
its sign. Z-DNA corresponds to a state of the molecule 
with the energy higher than the one in B form. Among the 
factors which induce B — Z transitions are {K — X)n — K 
peptide, polynuclear Pi-complexes, etc The inverse 
Z — B transition can be induced by KWGK peptide, 
daunomycin, and several other substances p^ . 

B — Z transitions have been extensively investigated 
that led to several models Recently rnolecular dy- 

namics and coarse-grained model [ll| simulations 
were carried out to provide with more satisfactory expla- 
nation of the molecular mechanisms of the B — Z transi- 
tions. In contrast, we here focus on the phenomenology 
of the transitions that are described by the minimalistic 
model presented above. The main goal was to demon- 
strate the possibility of transformations under the influ- 
ence of the external factors and temperature. 

As we mentioned it was shown that different helical con- 
formations of DNA may coexist as domains of the same 
long molecule Q. Now consider a chromatin, that is the 
complex of DNA and proteins in which the generic ma- 
terial is packaged inside the cells with nuclei [sj. Let us 
assume that a DNA molecule in chromatin happens to 



consist of right- and left-handed parts placed in a random 
order along the chain. Besides, chromatin structure is dy- 
namic jl7| that is important because changes in its struc- 
ture can be inherited independent of the DNA sequence 
itself, the so-called epigenetic inheritance. Our model is 
able, in principle, to mimick the case of DNA molecule 
with a dynamic structure as it provides with a plastic- 
ity mechanism via the fluctuating links (random bend- 
ing of DNA molecule caused by looping proteins) that 
(dis)connect base pairs. Thus a question arises if a DNA 
molecule in a chromatin can naturally exist in the state of 
matter we consider here, the pseudoscalar liquid crystal. 
This is truly an astonishing possibility and further ex- 
perimental and theoretical investigation of the molecular 
details would shed more light onto that. 

In conclusion, we have shown that a long DNA 
molecule, that consists of aperiodically ordered right- 
handed B and left-handed Z fragments, can form a novel 
condensed phase of matter, the pseudoscalar liquid crys- 
tal. We presented a minimalistic phenomenological model 
that describes the pseudoscalar liquid crystal phase of 
DNA and related phase transition phenomena. Gener- 
ally speaking, one can use the following ID Hamiltonian 
for theoretical treatment of the pseudoscalar liquid crystal 



dlh{l)x{l) 



with the continuous along the chain variable x(0 ^''^d. the 
interaction potential J{l,l'). Here x(0 pl^^ys the role of 
the liquid crystal chiral variable defined as n • V x n with 
n being the director for conventional liquid crystals 0], 
and h{l) being a heterogenous (random) external field. 

Further experimental research is needed to identify if 
this condensed state of matter exists in nature. Yet an- 
other foreseeable way towards the new state of matter 
would be assembling pseudoscalar liquid nano-crystals via 
DNA-programming 

We thank A. E. AUahverdyan, D. Mukamel, R.A. Roe- 
mer, E. I. Shakhnovich, and M. C. Williams for comments 
and discussions. This work was supported by National 
Science Council in Taiwan under Grant Nos. NSC 96- 
2911-M 001-003-MY3, NSC 96-2811-M-001-018 and NSC 
97-2811-M-001-055, by National Center for Theoretical 
Sciences in Taiwan, and by Academia Sinica (Taiwan) un- 
der Grant No. AS-95-TP-A07. 



[1] E. Schrodinger, What is life: The physical aspect of living 

cell (University Press, Cambridge, 1944). 
[2] J.D. Watson and F.H.C. Crick, Nature 171, 737 (1953). 
[3] B. Alberts et al, Molecular Biology of the Cell (4th edn. 

Garland Science, New York, 2002). 
[4] A. Wang et al, Nature 282, 680 (1979); A. Rich and S. 

Zhang, Nature Rev. Genet. 4, 566 (2003); P.C. Champ, S. 

Maurice, J.M. Vargason, T. Camp, and P.S. Ho, Nucleic 

Acids Res. 32, 6501 (2004). 



[5] S.C. Ha, K. Lowenhaupt, A. Rich, Y.G. Kim, and K.K. 

Kim, Nature (London) 437, 1183 (2005). 
[6] A.Yu. Grosberg and A.R. Khokhlov, Statistical Physics of 

Macromolecules, Chap. 3, §22 (AIP Press, Woodbury NY, 

1994). 

[7] Ya.B. Zel'dovich, Zh. Eksp. Teor. Fiz. 67, 2357 (1974) 

[Sov. Phys. JETP 40, 1170 (1975)]. 
[8] A pseudoscalar is a quantity that behaves like a scalar, 

except that it changes sign under a parity inversion. The 



examples of pseudoscalars are magnetic flux, spin helicity 

and the chirality which is considered in the present paper. 
[9] S. Chandrasckhar, Liquid Crystals (Cambridge University 
Press, New York, 1977); P.G. de Gennes and J. Frost, 
The Physics of Liquid Crystals (Clarendon Press, Oxford, 
1993). 

[10] L. Saiz and J.M.G. Vilar, Curr. Opin. Struct. Biol. 16, 344 
(2006); G.A. Maston, S.K. Evans, and M.R. Green, Annu. 
Rev. Genomics Hum. Genet. 7, 29 (2006); H.G. Garcia, P. 
Grayson, L. Han, M. Inamdar, J. Kondev, P.C. Nelson, R. 
Phillips, J. Widom, and P.A. Wiggins, Biopolymers 85, 
115 (2007). 

[11] J.V. Selinger and R.L.B. Selinger, Phys. Rev. Lett. 76, 58 
(1996); J.V. Selinger, Phys. Rev. Lett. 90, 165501 (2003). 

[12] A.E. Allahverdyan and K.G. Petrosyan, Europhys. Lett. 
75, 908 (2006). 



[13] D.A. Lavis and G.M. Bell, Statistical mechanics of lattice 

systems, v. 1 (Springer- Verlag, Berlin, 1999). 
[14] M.A. Fuortos, V. Cepeda, C. Alonso, and J.M. Perez, 

Chom. Rev. 106, 2045 (2006). 
[15] M.A. Kastenholz, T.U. Schwartz, and P.H. Hiinenberger, 

Biophys. J. 91, 2976 (2006). 
[16] W. Lim, Phys. Rev. E 75, 031918 (2007). 
[17] G. Felsenfeld and M. Groudine, Nature 421, 448 (2003). 
[18] D. Nykypanchuk, M.M. Mayo, D. van dor Lelic, and O. 

Gang, Nature 451, 549 (2008); S.Y. Park, A.K.R. Lytton- 

Jean, B. Lee, S. Wcigand, G.C. Schatz, and C.A. Mirkin, 

Nature 451, 553 (2008). 
[19] F.A. Aldaye, A.L. Palmer, and H.F. Sleiman, Science 321, 

1795 (2008). 



